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ABSTRACT 


The liquid and gas phase equilibria H 2 O + D 2 O = 2HD0 have 
been reinvestigated. By combining spectroscopic and vapor 
pressure data, fairly good agreement is found between the experi- 
mental and calculated enthalpy for the liquid phase reaction. 

It is found that the value of the equilibrium constant for the 
liquid phase reaction is very close to the geometric mean value 
of 4.0. 


THERMODYNAMIC PROPERTIES 
OF LIQUID H 2 0 AND D 2 0 AND THEIR MIXTURES 

By Mark Salomon 
Electronics Research Center 


SUMMARY 

The thermodynamics of the two isotopic equilibria 


H 2° (£) + D 2° (£) _ 2HD0 U) ^l] 

and 

H 2°(g) +D 2°(g) = 2HD0 (g> [2] 


have been studied with the object of obtaining more accurate 
values for the equilibrium constants. By using the experimental 
data for the vapor pressures of liquid H 2 O, D 2 O, and HDO, together 
with the spectroscopic gaseous data, the free energies and en- 
thalpies of reactions [l] and [2] were calculated. The only 
accurate experimental data relating to either [l] or [2] is the 
heat of mixing of liquid H 2 O and D 2 O. The results of the present 
calculation give a value of 33.2 calories per mole for the heat 
of mixing which is close to the experimental value of 32 calories 
per mole. 


INTRODUCTION 

Many investigators involved in hydrogen-deuterium studies 
in solution either directly or indirectly make use of the equi- 
libria 


H 2° ( l ) +D 2°U) ~ 2HD0 U) 


[ 1 ] 


Values for the equilibrium constant K(I) for reaction [l] have 
been measured directly (refs. 1-3) , have been calculated by an 
indirect method (ref. 4) , and have been calculated from spectro- 
scopic data (refs. 1,5-8). Despite these many attempts to arrive 
at a value of K(I) , there exist several discrepancies between 
calculated values , depending upon the source of the spectroscopic 
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data (refs. 9,10), the method of calculation, and also between 
the experimental values and the geometric mean value of four 
obtained by Pyper and Long (ref. 4) . In general, calculated 
values of K(I) based on Darling and Dennison's (DD) data (ref. 9) 
are higher than those based on the data of Benedict et. al. 

(ref. 10) (BGP) . The experimental values tend to indicate higher 
values than those obtained from the BGP data, but they are some- 
what smaller than those values obtained from the DD data. Weston 
(ref. 11) has recently discussed this problem and suggests that 
some accurate determination of K(II), 


HD 


(g) 


H o 0, . 
2 (g) 


= H 


2 (g) 


+ HDO 


(g) 


[ 2 ] 


would help resolve this problem. In this report, an additional 
criterion is suggested and investigated with the object of help- 
ing to resolve this problem. The additional criterion suggested 
here is the standard enthalpy of reaction [l], viz., AH°(lj which 
is chosen as a valid criterion because it is the only quantity 
which can be easily and accurately measured experimentally. The 
measured value of AH°(I) is 32 cal/mole (refs. 12,13). 


METHOD 

Since spectroscopic data for liquid species are not accurate 
enough to permit one to calculate K(I), the procedure used here, 
following Kirshenbaum ' s suggestion (ref. 7), is to relate reac- 
tion [l] to the gas phase equilibrium reaction: 


H 2°(g) + D 2°(g) = 2HD0 (g) 
via the relation 


K(I) 


K (III) 


p h 2 o * p d 2 o 
( P HDo) 


[3] 


[4] 


In Eq. [4], p° is the vapor pressure of the pure liquid indicated 
by the subscript. In terms of standard free energies, Eq. [4] 


2 



can be written as 


AG° (I) = AG° (III) - RT<2 In 


H 2° 

3 dho. 


- In 



[5] 


Because HDO can exist only in a mixture of 1^0 and D 2 O, its vapor 
pressure cannot be measured directly. This problem has often 
been avoided by assuming that the rule of the geometric mean 
(ref. 14) is applicable; i.e.: 


( p HDo) ( P H 2 o) ( P D 2 o) 

and hence that 

AG° (I) = AG° (III) . 


However, application of the rule of the geometric mean to the 
vapor pressures does not give acceptable AH 0 (I) values. The 
vapor pressure ratio 

__ 0 

Pli 2° 

PfiDO 


can be measured in dilute D-containing solutions by single-stage 
equilibrium experiments. Narten (ref. 8) has combined two sets 
of data (refs. 15,16) to arrive at the relation: 


^ O 

Pli 2° 

P HD0 


-0.033198 


11.6763 
T - 185.586 


[ 6 ] 


where T is temperature in degrees Kelvin. The experimental rela- 
tion for 


H„0 


P6 2 o 
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is given by Miles and Menzies (ref. 17) as 


log 


H„0 


D 2° 


16.9987 - — 8 ^, 843 - 7.4972 log T 
+ 9.7611xlO -3 T 4.4288 x10 -6 T 2 . 


[ 7 ] 


Another relationship obtained by Whalley (ref. 18) is that result- 
ing from a least squares refinement of several sets of data and is 


log 


P «2° 


-0.030661 + 


9.14056 
75.753 + t 


[ 8 ] 


where t is the temperature in degrees centigrade. Narten (ref. 8) 
used Eqs . [6] and [8] to evaluate the term in curly brackets in 

Eq. [5] and combined this with the free energy functions cal- 
culated by Friedman and Haar (ref. 19) and found AH° (I) = 137 
cal/mole. Whalley (ref. 18) has pointed out that Eq. [8] gives 
acceptable results for temperatures above 100°C, but deviates 
significantly from the experimental relation [7] below this 
temperature. In the present calculations , the correct value of 
AH 0 (I) was obtained only by using Eqs. [6] and [7]. 

To calculate K(I) values, the corresponding K(III) values 
are required. To calculate K(III) values, the spectroscopic 
data of BGP and DD [as quoted by Kirshenbaum (ref. 7)] were used. 
The familiar statistical mechanical procedures (ref. 20) were 
employed in the calculation. The translational partition func- 
tion (p.f.) was calculated in the usual manner. Since there is 
no explicit formula for the rotational levels of an asymmetric 
top molecule, use was made of the approximate relation given by 
Herzberg (ref. 21): 


^r 


exp 


(BC) 1 ^ 2 hc 
4kT 



[9] 


where A, B, and C are the effective rotational constants defined 
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by 


A = A e - 2>i (v.+|) 

B = B e - ( v i + |) 

C = C e - (v. + ±) . [10] 

BGP have evaluated the nine molecular constants ai for each iso- 
tope whereas DD have evaluated similar, but not equivalent, a-j_ 
values. Herzberg (ref. 21) re-evaluated the DD aj_ values in 
terms of the constants for H 2 O only, and they are very simi- 
lar to those reported by BGP. Hence the BGP ct± values have been 
used for HDO and D 2 O in all the calculations presented in this 
report . 

The vibrational p.f. was evaluated from the relation 

q v = 5 exp - G(v ± )|^ [11] 


where 


G(v i ) 


£ 


w i ( v ± + 1 ) 


+ £ 




[ 12 ] 


In Eq. 12, the coj_ are the zero-order frequencies and the xj_j are 
the anharmonicities . Both the u)i and xij values reported by BGP 
and DD differ slightly except for HDO where the difference is 
very large. The result is that the zero point energies for HDO 
vary accordingly as seen in Table I. The DD data indicate the 
only important states to be (0, 0, 0) and (0, 1, 0) and the BGP 
data show no important contributions to G(v^) beyond the second 
vibrational state. The calculations presented here, however, 
were carried out for all possible combinations of V]_, V 2 , and V 3 
up to the fourth vibrational state. The results were practically 
identical whether the first, second, third, or fourth vibrational 
states were considered. All these calculations were carried out 
on a Univac 1108 computer located at the NASA Manned Spacecraft 
Center, Houston, Texas. In all cases the effects due to nuclear 
spin have been ignored. 
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TABLE I 


ZERO POINT ENERGIES* FOR ISOTOPIC WATERS 



h 2 o 

hdo d 2 o 

Data of Dennison 
Data of Benedict 

4631.247 

4634.253 

4018.585 3386.238 

4032.162 3388.687 


*In cm 


RESULTS AND DISCUSSION 

Table II gives the results of the calculations for the 
gaseous reaction [3]. Table III gives the data for the liquid 
phase reaction [l] obtained from the vapor pressure relations 
[6] and [7] and Eq. [5]. Table IV gives the results of a least 
squares fit to the relation 


AG° (I) 


A + BT + CT 


2 


The first entry in Table IV, which is one of the ERC results, is 
in very good agreement with the experimental AH 0 (I) value. 
Kirshenbaum 1 s K(III) values also led to an acceptable value for 
the enthalpy. All calculations using the BGP data did not give 
very good AH° (I) values. It should be noted that the direct 
determination of K(I) by Friedman and Shiner (ref. 2) gives a 
very good value for AH° (I) . 


TABLE II 


THERMODYNAMIC QUANTITIES FOR REACTION [3] 


°C 

Data of Dennison 

Data of Benedict 

K (III) 

AH° (III) * 

K (III) 

AH° (III)* 

0 

25 

50 

75 

80 

100 

3.985 

4.019 

4.047 

4.071 

4.075 

4.083 

56.62 

54.61 

52.41 

50.05 

49.56 

48.561 

3.409 

3.471 

3.524 

3.568 

3.576 

3.592 

118.10 
116 . 06 
113.82 
111.41 
110.90 
109 . 89 


*Units of AH are calories per mole. 
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TABLE III 


THERMODYNAMIC QUANTITIES FOR REACTION [l] 



Data of 

Dennison 

Data of 

Benedict 

c 

K(I) 

> 

E 

O 

H 

* 

K ( I ) 

AH° (I) 

0 

4.010 

15.33 

3.432 

59.82 

25 

4.010 

33.15 

3.463 

94.54 

50 

4.035 

67.76 

3.522 

129.19 

75 

4.086 

! 106.51 

3.582 

167.98 

80 

4.094 

114.95 

3.593 

176.43 

100 

4.108 

132.60 

3.614 

194.10 


*Units of AH are calories per mole. 


TABLE IV 


POWER SERIES CONSTANTS A, B, C FOR REACTION [l] 


AH 0 (I) * at 
298 . 15°K 

A 

B 

C 

Notes/ 

References 

33.15 

-174.32 

-1.486 

-2.322 

10" 3 

a 

94.54 

-113.15 

-1.399 

-2.324 

10~ 3 

b 

39 . 27 

-173.61 

-2 . 384 

-6.370 

10~ 4 

2 (g) 

21.06 

-9.56 

-2.613 

-3.146 

10" 4 

6 

36.09 

-220.43 

-1.052 

-2 . 877 

10~ 4 

7 

65.81 

-358.00 

-2.506 

-4.101 

10' 3 

2f 

84 . 94 

-98.79 

-1.546 

-2.067 

10~ 3 

8 

100.63 

230.71 

-3.977 

1.459 

10~ 3 

c 

162.09 

292.47 

-3 . 894 

1.462 

10~ 3 

d 

50.07 

73.24 

-3.051 

2.387 

10- 4 

6 (e) 


*AH units are calories per mole. 


a, c = results of present calculations using the DD data 

b, d = results of present calculations using the BGP data 

c, d, e = results based on Eqs . [5], [6], and [8] 

g, f = ref. 2; experimental and theoretical data, 

respective ly 

All other entries were calculated from Eqs. [5] - [7], 
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The results of the present calculations can be extended to 
additional isotopic equilibria. To accomplish this, the p.f. 
ratio 


q HDO 

qH 2° 


for the liquid species is needed. This ratio can be calculated 
from the relation: 


q HDO 
q H 2 0 a 



[13] 


where the ratio 


q H 2 0 £ 


has been calculated previously in the temperature range of 
interest (ref. 22) . 

Two important equilibria are 


2° (g ) 

+ HD, . 

(g) 

■ HD0 (g) 

+ H 2(g) 

2° (£) 

+ HD, . 

(g) 

- HD V) 

+ H 2 (g) 


The p.f. ratio 



[14] 

[15] 


was obtained as a function of temperature from an NBS monograph 
(ref. 23) , whereas the ratios for the water species were obtained 
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from the present calculations resulting from the application of 
the DD data. K(14) and K(15) values are given in Table V as a 
function of temperature. 


TABLE V 

EQUILIBRIUM CONSTANTS K(14) AND K(15) 


°c 

0 

25 

50 

75 

80 

K(14) 

4.340 

3.752 

3.360 

2.986 

2.929 

K (15 ) 

4.780 

4.000 

3.449 

3.049 

3.982 


CONCLUSIONS 

The present work was performed with the intent of clarify- 
ing some of the existing uncertainties in the value of K(I) and 
the problems associated with its calculation. The ERC calcula- 
tions agree very well with the measured quantity AH°(I) and 
with Pyper and Long's value at 25°C. The ERC results are also 
in fair agreement with the measurements of Friedman and Shiner 
(ref. 2) and hence lend support to the validity of the higher 
value of K (I ) . 

Despite the fact that the BGP study (ref. 10) was much more 
extensive than the earlier one by Darling and Dennison (ref. 9) , 
it gives very low K(I) values and very high AH°(I) values. The 
reason for this difference has been traced (refs. 4, 11) to 
differences in anharmonicity constants for HDO which results in 
significant differences between the zero-point energies for the 
gaseous molecules. Since the HDO spectrum, by necessity, con- 
tains equilibrium quantities of H 2 O and D 2 O, it is quite under- 
standable that the frequencies and molecular constants for HDO 
are the least reliable (cf . discussions in references 4 and 11) . 

Darling and Dennison (ref. 9) did not observe the HDO spec- 
trum but calculated the required molecular constants by the 
usual methods based on the H 2 O and D 2 O spectra. The direct 
summation of the vibration energies does not account for the 
present success in the calculated value for AH° (I) since this 
method has also been applied to the BGP data. 

Weston (ref. 11) discussed the use of other equilibria which 
would aid in the choice between the more accurate spectroscopic 
data. For example, he cites the isotopic equilibria [14] for 
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which the foregoing equilibrium constant has been calculated 
here. Weston gives two references (refs. 24,25) in which K(14) 
was indirectly estimated to be 3.41 and 3.54 at 25°C. These 
values are in good agreement with those values calculated using 
the BGP data (ref. 11) but are not in particularly good agree- 
ment with the ERC value of 3.75 based on the DD data. According 
to Kirshenbaum (ref. 7) , however, Grosse (ref. 26) measured a 
value of 3.88 for K(14) at 25°C. Thus this important point needs 
further experimental confirmation. 

While the present calculations give acceptable results, it 
should be pointed out that Wolfsberg (ref. 27) has questioned 
the use of Eq. [ 1 l] , since he questions the usual assumption 
that the anharmonicity constants are invariant to isotopic sub- 
stitution. If this is true, Eq. [ 11 ] cannot be used to calculate 
the zero-point energies. 


Electronics Research Center 

National Aeronautics and Space Administration 
Cambridge, Massachusetts, January 1969 
120-34-02-01 
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